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TIME  CONSTRAINTS  ON  MEASURING 
BUILDING  R-VALUES 


Stephen  N.  Flanders 


INTRODUCTION 

This  report  describes  the  time  constraints  on  the 
accuracy  of  measuring  thermal  resistance  (R-values)  of 
building  components.  It  gives  a  method  for  estimating 
the  characteristic  time  constant  of  building  elements 
and  demonstrates  the  magnitude  of  error  that  different 
weather  events  may  introduce  into  thermal  measure¬ 
ments.  This  report  should  help  an  investigator  to  de¬ 
termine  the  appropriate  duration  of  measurement  for 
determining  R-values  and  demonstrate  the  magnitude 
of  error  in  the  resulting  data. 

A  high  R-value  helps  to  slow  down  heat  flow  out  of 
a  building  and  thereby  conserves  space  heating  energy. 
Unfortunately,  the  thermal  performance  of  a  building 
is  usually  worse  than  the  design  thanks  to  incorrect  in¬ 
stallation  and  deterioration  of  insulation.  Therefore, 
if  we  want  to  determine  the  actual  R-value  of  building 
components,  we  must  either  remove  samples  to  a 
laboratory  or  measure  the  components’  thermal  properties 
in  place. 

A  laboratory  test  can  accurately  determine  building 
component  thermal  properties  but  sample  removal 
damages  the  building  and  does  not  reflect  on-site  con¬ 
ditions.  In  the  laboratory  we  can  maintain  constant 
temperature  and  heat  flow  and  determine  R-values  ac¬ 
cording  to  a  basic  heat  flow  equation: 

R  =  LT/q  (1)  ' 

where 

at*  rw-rc 

R  -  thermal  resistance 
Tw  =  temperature  on  warm  side  of  sample 


rc  =  temperature  on  cold  side  of  sample 
q  =  heat  flux. 

As  long  as  the  sample’s  thermal  conditions  maintain  a 
steady  state,  sensors  in  the  test  apparatus  give  meaning¬ 
ful  values  for  the  two  temperatures  and  the  rate  of  heat 
flow. 

Field  measurement  of  R-values  doesn’t  damage  the 
building,  but  it  doesn’t  give  us  steady-state  conditions 
either.  If  we  measure  the  R-value  of  a  wall,  we  apply 
temperature  and  heat  flow  sensors  on  the  indoor  surface 
and  a  temperature  sensor  at  a  corresponding  location  on 
the  outdoor  surface.  Figure  1  shows  a  strip  chart  re¬ 
corder  monitoring  the  output  of  such  sensors  for  several 
locations  on  a  wall.  Fluctuating  temperatures  on  either 
side  of  the  wall  prevent  it  from  passing  heat  flow  at  a 
steady  rate,  so  that  at  any  one  time  the  measured  tem¬ 
perature  difference  across  the  wall  (A T)  and  heat  flux 
(q)  do  not  render  accurate  R-values  in  eq  1 .  Therefore, 
a  major  problem  of  on-site  measurement  is  to  determine 
the  minimum  duration  of  measurement  to  obtain  an 
accurate  R-value. 

This  problem  occurs  because  a  building  component 
not  only  resists  the  flow  of  heat,  but  also  stores  heat. 
Each  building  material  has  its  own  resistance  to  heat 
flow  (R-value),  heat  storage  capacity  (specific  heat)  and 
density.  These  factors  contribute  to  the  time  for  a 
change  in  temperature  across  a  building  component  to 
reach  a  new  steady-state  value.  If  AT  across  the  com¬ 
ponent  increases,  the  component  stores  heat  up  to  its 
capacity  at  that  temperature  potential  until  it  passes 
as  much  through  the  cool  side  as  enters  througi  the 
warm  side.  If  the  temperature  difference  decreases,  the 
component  gives  off  heat  until  it  reaches  its  new  dimi¬ 
nished  capacity  and  again  passes  as  much  as  enters.  The 


Figure  1.  On-site  measurement  of  building  R -values. 
Sensors  under  the  paper  on  the  wall  measure  heat  flow 
and  temperature  on  the  Indoor  surface  while  a  temper¬ 
ature  sensor  outdoors  completes  the  needed  Information 
which  a  strip  chart  records  for  several  locations. 

length  of  time  required  for  the  component  to  reach  a 
steady-state  value  depends  on  its  characteristic  time 
constant. 

With  these  basic  principles  in  mind,  we’ll  see  that 
the  accuracy  of  field  measurements  of  building  R-values 
depends  on  three  main  factors: 

1 .  Temperature  difference  (AT)  across  the  component. 

2.  Variation  in  AT  prior  to  and  during  the  observa¬ 
tion  period. 

3.  Characteristic  time  constant  of  the  component. 
Sensor  accuracy  and  the  impact  of  measurement  on 
altering  the  actual  conditions  are  the  subjects  of  a  col¬ 
league’s  unfinished  study. 

The  larger  the  average  AT  across  a  component,  the 
less  a  given  variation  in  AT  will  affect  the  measurement. 
We  can  classify  each  of  die  many  possible  causes  for 
variation  in  AT  as  either  cyclical  or  random.  The  pas¬ 
sage  of  seasons  and  days  causes  cyclical  temperature 
changes.  However,  seasons  change  so  gradually  that 
they  are  essentially  constant  for  the  time  we  would 
measure  an  R-value.  We  can  account  for  the  values  of 
AT  that  cause  the  heat  to  flow  and  all  the  values  of 
q  that  represent  the  amount  of  heat  flow  in  a  new 


formula  for  R-value: 

n  In 

r = £  We  «■  w 

s=i  /  i=i 

where  AT;  is  the  difference  in  temperature  between  the 
inside  and  outside  at  the  i-th  measurement  and  q,  the 
heat  flow  at  the  same  measurement  time.  This  equa¬ 
tion  works  well  for  cyclical  variations  in  temperature, 
as  long  as  we  measure  frequently  for  one  full  cycle,  but 
has  a  drawback  for  random  changes. 

If  random  temperature  changes  are  significant  before 
or  during  the  time  we  measure,  we  then  must  be  more 
careful  with  our  use  of  eq  2.  We  should  measure  for 
a  long  enough  time  period  to  know  that  the  heat  taken 
into  or  rejected  from  storage  in  the  building  is  not  large 
compared  to  the  amount  we  saw  pass  through.  The 
component’s  thermal  time  constant  determines  how 
much  of  this  heat  is  unaccounted  for. 

A  component’s  time  constant  tc  comes  from  the 
thermal  properties  of  resistivity,  specific  heat,  density 
and  thickness  for  each  type  of  building  material  in  a 
wall  or  roof  layer.  The  value  tc  helps  us  to  understand 
how  quickly  the  component  responds  to  changes  in 
AT-. 

<3) 

where 

re  =  characteristic  time  constant  of  building  com¬ 
ponent 

<?n  =  (a>k)*  ,  conversion  constant  adjusting 
thickness  of  layer  to  make  material  uniform 
throughout  wall 

a„  =  rnC„dn,  reciprocal  of  diffusivity  of  n-th 
layer 

r„  =  resistivity  of  n-th  layer  (published  value) 

Cn  =  specific  heat  of  n-th  layer  (published  value) 
dn  -  density  of  n-th  layer  (published  value) 
xn  =  thickness  of  n-th  layer 
«l(  =  a„  at  layer  k  chosen  for  normalization. 

This  equation,  derived  in  Appendix  G,  transforms  the 
thermal  properties  of  each  layer  into  those  of  one  com¬ 
posite  material,  but  adjusts  the  thickness  of  each  layer 
to  compensate  and  retain  each  layer's  original  time  be¬ 
havior  so  that  we  may  calculate  fc  as  if  it  represented 
a  single  material. 

With  tc  we  can  predict  how  rapidly  a  building  element 
responds  to  changes  in  AT  because  te  reflects  both 


2 


the  resistance  to  heat  flow  and  the  heat  storage  capacity 
of  the  component.  If  A7"  changes  linearly  (like  a  ramp), 
then  q  will  follow,  delayed  by  a  period  tc.  If  A7"  changes 
abruptly  (like  a  step),  then  q  follows  with  time  t  ac¬ 
cording  to  and  reaches  63.2%  of  its  new  equil¬ 

ibrium  value  at  t  =  tc.  In  the  next  section,  well  see 
that  this  value  comes  from  an  expression  of  Carslaw 
and  Jaeger  (1959).  The  time  constant  also  controls 
the  component’s  response  to  cyclical  changes. 

Such  delays  in  thermal  response  occur  because  the 
component  temporarily  gives  off  or  stores  more  heat 
than  it  can  over  the  long  term;  tc  is  the  key  to  the  mag¬ 
nitude  of  error  that  this  transient  phenomenon  causes. 
Equation  3  may  appear  to  be  academic  because  the 
variables  are  experimental  unknowns.  Nevertheless, 
inspection  of  the  component  or  of  construction  drawings 
together  with  published  data  about  thermal  and  material 
properties  will  give  sufficient  accuracy  for  field  mea¬ 
surement. 

We  can  get  a  feel  for  some  typical  building  wall  time 
constants  tc  using  eq  3  and  data  from  ASHRAE  (1977) 
in  Appendix  B.  A  2x4  insulated  frame  wall  has  a  tc 
of  1.25  hr,  a  2x6  wall  1.78  hr  and  and  an  insulated  4-in. 
brick  and  8-in.  concrete  block  masonary  wall  2.82  hr. 
Later  well  use  these  values  in  some  examples,  but 
first  well  discuss  the  contributions  of  other  investi¬ 
gators  to  the  problem  of  determining  the  time  con¬ 
straints  on  thermal  measurements. 


FIELD  MEASUREMENT  AND  ANALYSIS  OF 
TRANSIENT  HEAT  FLOW 

Much  literature  on  the  measurement  of  heat  flow 
under  uncontrolled  conditions  comes  from  the  study 
of  soil  temperatures.  For  example,  Lachenbruch 
(1959)  discusses  how  to  analyze  the  seasonal  variation 
of  temperature  in  stratified  soils. 

Relatively  little  is  written,  however,  about  how  long 
to  use  heat  flux  sensors  and  thermocouples  in  measuring 
the  thermal  performance  of  building  components,  al¬ 
though  this  is  a  common  investigation.  Lorentzen  et  al. 
(undated)  outline  heat  flux  and  temperature  sensor 
application  methods  and  limitations.  Peavy  et  al. 

(1975)  and  Burch  (1976)  of  the  National  Bureau  of 
Standards  are  two  frequent  users  of  this  approach  to 
in-situ  measurement  of  building  components.  Such 
research  often  recognizes  time  constraints  by  measuring 
for  ample  durations  or  24-hour  periods. 

Poppendieck  et  al.  (1976)  outline  some  basic  con¬ 
siderations  in  non-steady-state  measurement  of  heat 
transfer  in  buildings.  They  give  the  mathematical  solu¬ 
tions  for  temperature  as  a  function  of  depth  and  time 
in  an  idealized  one-material  wall  for  both  a  sinusoidal 


and  a  step  change  in  temperature  across  the  wail.  They 
calculate  the  phase  lag  from  a  sinusoidal  variation  in 
temperature  across  a  2-in.-thick  cork  board  to  be  1 1 
hours  and  then  demonstrate  that  measurement  of  periodic 
temperature  difference  and  heat  flow  for  one  period 
should  result  in  a  ratio  equal  to  the  steady-state  R-value. 

A  single  laboratory  test  using  this  2-in.  cork  panel  in 
a  guarded  hot  box  corroborated  their  theoretical  ex¬ 
pectations  for  R-value  and  time  constants. 

Poppendieck  et  al.  (1976)  suggest  that  measurement 
over  a  full  24-hour  cycle  will  give  accurate  data  for  the 
R-value.  But  they  do  not  answer  how  to  cope  with  an 
underlying  trend  for  changing  temperature  or  for  a 
change  of  amplitude  in  temperature  extremes.  They 
also  suggest  that  a  nighttime  reading  can  suffice  for 
walls  less  than  R  =  10.  However,  since  specific  heat, 
density  and  wall  thickness  also  influence  the  length  of 
time  for  a  building  component  to  stabilize  its  thermal 
behavior,  this  recommendation  lacks  authority. 

Further  consideration  of  the  time  constraints  on 
thermal  measurement  has  come  from  1 )  mathematical 
analyses,  2)  electrical  and  hydraulic  analogs  or  3)  com¬ 
puter  models. 

Carslaw  and  Jaeger  (1959)  give  a  mathematical  solu¬ 
tion  for  a  slab  equilibrating  to  constant  temperatures 
on  each  side.  Their  equation  demonstrates  the  basic 
qualities  of  a  component  layer: 


where 

v  =  temperature  at  x 

<7n  =  coefficient  in  a  Fourier  series  representing 
initial  conditions 
K  =  thermal  diffusivity 
/  =  thickness  of  slab 
x  =  distance  into  slab 
t  =  time. 

The  sine  term  pertains  to  the  location  in  the  slab.  The 
exponent  is  the  basis  for  eq  3  and  describes  the  time- 
related  behavior  of  the  slab.  Evaluating  v  to  n  =  1  is 
usually  sufficiently  accurate. 

Granholm  (1971),  Trethowen  (1972)  and  Petzold  et 
al.  (1974)  discuss  mathematical  techniques  for  analyzing 
non-steady-state  heat  flow  through  walls.  Ullah  et  al. 
(1976)  refine  the  multiple  harmonic  Fourier  method  to 
calculate  periodic  heat  flow  in  building  walls.  Sonderegger 
(1977)  analyzes  the  relative  advantage  of  having  a  mas¬ 
sive  outside  layer  with  insulation  inside  the  building 
envelope,  or  vice-versa,  with  thermal  response  factors 
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based  on  Fourier  transforms.  Peavy  (1 978)  develops 
a  streamlined  method  for  calculating  thermal  response 
factors  that  includes  the  effects  of  convection  and  radia¬ 
tion  where  only  conduction  is  usually  considered.  Un¬ 
fortunately,  such  mathematical  analyses  cannot  readily 
handle  typical  fluctuations  in  A T  across  the  building 
walls. 

Electrical  and  hydraulic  analogs  have  also  been  used 
to  model  the  thermal  properties  of  building  walls.  Hawk 
and  Lamb  (1963)  come  closest  to  representing  the  prob¬ 
lem  of  describing  the  transient  thermal  behavior  of 
building  walls  with  a  hydraulic  model  employing  reser¬ 
voirs  and  valves  to  represent  heat  storage  capacity  and 
thermal  conductivity  of  building  layers. 

Ultimately,  computer  models  are  the  most  versatile 
of  these  methods  for  gaining  insight  for  accommodating 
unsteady-state  conditions  in  field  measurements.  The 
finite  difference  method  (Forsythe  1960)  for  approxi¬ 
mating  differential  equations  describing  the  heat  balance 
for  selected  locations  in  ihe  wall  is  quite  common. 

Fourier  transforms  discussed  above  and  finite  element 
techniques  (Aziz  1 972)  are  also  current. 

A  CLOSER  LOOK  AT  HANDLING  THE  CONSTRAINTS 

Measurements  spanning  a  long  duration  will  mini¬ 
mize  the  effect  of  transient  errors  on  R-value$.  How¬ 
ever,  anyone  taking  a  measurement  wants  to  spend  the 
least  time  possible  and  yet  know  that  the  error  is  ac¬ 
ceptable.  The  magnitude  of  error  is  time-dependent  on 
the  nature  of  the  building  element  and  the  climatic 
changes  that  occur  during  measurement. 

Let’s  think  of  the  difference  between  inside  and 
outside  temperatures  A T  as  corresponding  to  a  potential 
heat  flux.  Actual  heat  flux  will  tend  towards  the  value 
of  this  potential  until  AT  changes.  The  true  potential 
heat  flux  value  is  qa  =  A7//?a,  where  /?a  is  the  actual 
R-value  of  the  wall.  In  practice,  we  assume  a  value  of 
R,  based  on  theoretical  calculation  or  prior  measure¬ 
ment.  This  gives  heat  flux  potential  as  qp  =  A T/R. 

This  is  a  convenient  way  to  convert  temperature  readings 
into  units  of  heat  flow  and  compare  the  lagged  behavior 
of  the  actual  heat  flow  being  measured. 

Next,  we'll  compare  the  difference  between  the  area 
under  curves  representing  the  potential  flux  qp  and  the 
measured,  delayed  response  in  heat  flux  qr  over  the 
measurement  period  rb  and  determine  its  significance 
as  an  error  of  measurement.  In  each  case  well  find  the 
error  percentage  from 


Cyclical  change 

Development  of  eq  2  demonstrates  why  measure¬ 
ment  over  one  full  period  (such  as  1  day)  will  render 
a  good  calculation  of  R  for  true  cyclical  change.  Let’s 
assume  that  AT  fluctuates  according  to  a  cyclical  curve, 
such  as  the  one  at  the  top  of  Figure  2,  and  that  q  fluc¬ 
tuates  according  to  any  other  curve  of  the  same  period 
length  but  out  of  phase  (as  the  figure  shows).  Inspection 
of  the  areas  enclosed  by  Ihe  two  curves  for  any  full 
cycle  shows  each  area  to  be  constant.  Equation  2,  trans¬ 
lated  into  integrals,  is  the  ratio  of  these  areas  and  renders 
R  accurately.  Any  measurement  that  does  not  represent 
a  full  cycle  is  in  error  to  the  extent  that  the  ratio  of 
the  areas  is  disproportionate. 

The  diurnal  cycle  is  not  strictly  sinusoidal.  However, 
an  examination  of  the  properties  of  a  sinusoidal  signal 
passing  through  a  building  component  will  give  some 
insight  into  the  nature  of  the  response  curve.  The  de¬ 
laying  effect  of  a  building  component  may  alter  the 
heat  flux  response  in  two  ways:  1 )  the  measured  heat 
flux  qt  may  lag  behind  the  potential  heat  flux  qp  by  a 
fixed  interval  but  with  the  same  period,  and  2)  qTlqp 
may  decrease.  This  is  illustrated  in  Figure  3. 

Phase  shift  and  attenuation  vary  with  the  time  con¬ 
stant  of  the  building  component  and  the  period  of  the 
signal.  One  other  factor  is  important:  whether  the  re¬ 
sistance  to  heat  flow  within  the  wall  is  relatively  great 
compared  to  the  resistance  into  and  out  of  the  wall. 

The  greater  the  relative  internal  resistance,  the  higher 
the  order  of  delay  it  represents  in  Figures  4  and  5. 

These  figures  from  Forrester  (1961 )  depict  the  relation¬ 
ship  between  the  time  constant  of  a  wall  and  the  at¬ 
tenuation  of  amplitude  of  heat  flux  with  the  corresponding 
lag  in  flux. 

Since  24  hours  is  the  most  likely  cycle  period  to 
influence  a  wall  measurement,  Table  1  represents  ror 
this  time  period  the  time  ratios,  amplitude  ratios  and 
response  lag  for  our  wall  examples  as  shown  on  the  cover. 

The  percentage  of  measurement  error  during  the  course 
of  sinusoidal  variation  in  AT,  given  arbitrary  starting 
and  ending  times  for  measurement,  comes  from 

£  (12/rr)  (aU-b  V) 

*  c(f2-f,  )-(12e/ir)  (U) 


(6) 


Figure  2.  Two  arbitrary  curves  of  the  same  periodicity  tp 
but  with  one  lagged  behind  the  other.  The  ratio  of  the 
shaded  areas  of  each  curve  will  always  be  constant,  as  long 
as  the  measurement  duration  equals  tp. 


Figure  3.  Sinusoidal  change  in  AT  and  corresponding  qp  (< ^original 
amplitude  of  signal,  b-attenuated  amplitude  of  response).  Mea¬ 
sured  qr  lags  behind  qp  by  tL  and  attenuates  by  a  ratio  of  b/a. 


where 

U  =  cos  (nt}  /I 2)-cos(jrf 2/1 2) 

V-  cos  (ff/12)  (f,-ft)-cos  (tt/1  2)  (ti-tL) 

=  error  percentage  for  sinusoidal  change 
a  -  amplitude  of  potential  heat  flow  variation 
from  mean 

b  -  amplitude  of  heat  flow  response  variation 
from  mean 

c  =  mean  value  of  heat  flow 
f|  =  time  measurement  begins 
=  time  measurement  ends. 

This  equation  is  derived  in  Appendix  D. 

A  few  applications  of  eq  6  can  quickly  convince  us 
to  adhere  to  a  measurement  duration  of  approximately 
one  full  period.  Calculations  plotted  in  Figure  6 


demonstrate  that  large  errors  in  measurement  can  oc¬ 
cur,  depending  on  the  phase  of  the  cycle  during  the 
course  of  measurement  and  when  the  measurement 
begins  and  ends. 

The  figure  shows  the  thermal  response  of  our  sample 
walls  to  a  sinusoidal  variation  in  AT.  The  percentage 
error  in  measurement  is  plotted  as  a  function  of  the 
duration  of  measurement  according  to  the  following 
conditions: 

1 .  Mean  outdoor  temperature  is  either  40°F  with 
a  swing  of  ±30°F  or  30°F±10°F. 

2.  Indoor  temperature  is  70°  F. 

3.  Measurement  begins  at  various  times  (f  ,  =  0,1 ,6, 
12,18  hr)  after  sin(0). 

A  sample  calculation  demonstrates  how  to  obtain  a 
value  on  the  graph  from  eq  6.  Assume  a  2x6  frame 
wall  with  R  =  20  hr  ft2  °F/BTU,  amplitude  ratio  =  0.93 
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Figure  4.  Ratio  of  output  amplitude  to  Input  amplitude  vs  ratio  of  delay  time  to  input  period.  * 
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Figure  5.  Phase  shift  to  exponential  filter  vs  ratio  of  delay  to  period.  * 


•Reproduced  with  permission  from  Forrester  (1961). 
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Figure  6.  Percentage  error  from  a  sinusoidal  change  In  AT  as  a  function  of  measurement  duration  for  measurement 
start  time  r,  *  0, 1,6, 12,1 8  hr.  The  left-hand  graphs  represent  an  average  temperature  of  40° F  with  a  variation  of  ±3(fF. 
The  left-hand  graphs  represent  3OUF±10°F.  The  top  pair  of  graphs  show  the  behavior  of  a  2x4  frame  wall  (R  =  14, 
amplitude  ratio  (h.r)  =  0.96,  tL  *  1.25)  the  middle  pair  a  2x6  frame  wall  (R  =  20,  Ar  *  0.93,  tL  =  1. 78)  and  the  lower 
an  Insulated  masonry  wall  (9.  =  1 1,  Ar  =  0.92,  tL  =  2. 77). 


and  response  lag  =  l  .78  hr  (from  Table  1).  What  will 
the  error  of  measurement  be  after  measuring  for  5  hours 
after  the  1 2-hour  point  in  the  40°±30°F  sinusoidal  cycle? 
From  q  =  AT/R,  a  =  (70-40+30-70+40+30, /2{20),  b 
-  0.93u  and  c  =  (70-40)/20.  Therefore  substituting 
U-  cos  (ir/12)(12)-cos  (jt/12)(12+5)  and  y  =  cos 
(ff/1 2)  (1 2-1 .78)-cos  (ir/1 2)(1 2+5-1 .78)  into  eq  6 
renders  Ey,  =  -25.8%  in  agreement  with  the  circled 
data  point  on  the  left  middle  graph.  Programs  for  the 
Hewlett  Packard  HP-25  hand  calculator  that  calculate 
E%  for  different  changes  in  A T  appear  in  Appendix  E. 

The  40°±30°F  case  represents  an  extreme.  The  mean 
A r  could  provide  adequate  measurement  accuracy  if 
it  were  steady-state,  but  the  variation  in  A7"  which 
might  occur  with  a  southern  exposure  on  a  clear  spring 
day  in  New  England  makes  the  duration  of  measure¬ 
ment  crucial  for  accuracy.  The  30°±1 0°F  case  offers 
a  more  favorable  mean  AT,  and  together  with  the  dimin¬ 
ished  variation,  makes  the  error  less  than  5%  for  ar¬ 
bitrary  measurement  periods  longer  than  1 5  hours. 

The  graphs  in  Figure  6  clearly  show  how  percentage 
error  increases  with  increasing  time  constants  of  building 
elements. 

Since  data  from  diurnal  cycles  are  not  sinusoidal,  they 
do  not  necessarily  adhere  closely  to  the  guideline  for 
error  determination  that  eq  6  represents.  However,  the 
peaks  and  valleys  of  the  input  AT  values  often  corre¬ 
spond  clearly  with  extremes  in  q,  by  a  fixed  lag  time. 
Integrating  the  data  for  both  variables  for  corresponding 
lagged  periods  should  diminish  the  error  of  measure¬ 
ment  for  less  than  one  full  cycle’s  monitoring. 

Random  change 

While  the  significant  cyclical  events  hinge  around  a 
diurnal  cycle,  random  changes  in  temperature  may  be 
sudden  or  gradual.  The  most  sudden  event  may  be  a 
cloud  passing  across  the  sun,  or  a  weather  front  may 
change  air  temperature  quite  abruptly.  Generally, 
however,  air  temperature  changes  gradually  for  a  period 
after  frontal  passage.  These  and  many  other  events  can 
affect  the  change  in  AT  across  a  building  component 
and  be  an  underlying  disturbance  to  the  accuracy  of 
measurement. 

Abrupt  change  in  outside  temperature  might  cause 
temperature  and  heat  flow  within  a  building  component 
to  respond  as  in  Figure  7.  An  immediate  reading  of 
AT  and  heat  flux  would  give  a  misleading  calculated 
R-value  from  eq  1 .  Therefore,  we  would  want  to  mea¬ 
sure  for  an  R-value  long  enough  that  the  amount  of 
heat  the  component  stores  or  releases  when  the  en¬ 
vironment  warms  or  cools  is  not  large  compared  with 
the  total  flow  through  the  wall. 

An  abrupt  (step)  change  or  a  linear  (ramp)  change 
in  AT  affects  the  accuracy  of  eq  1 .  Since  nature 


produces  a  combination  of  such  signals,  we  should  ap¬ 
preciate  the  relative  impact  of  each. 

Step  change 

A  cloud  passing  overhead  can  cause  a  rapid  change 
in  surface  temperature  on  a  wall  oriented  to  the  sun.  In 
less  than  an  hour  a  cold  front  can  cause  a  rapid  change 
in  temperature  that  persists  long  afterwards.  If  tc  of 
the  building  component  is  large  enough,  either  event 
may  approximate  an  abrupt  change  in  temperature  as 
depicted  in  Figure  8.  This  step  change  creates  a  poten¬ 
tial  flux  qp  that  the  actual  flux  q,  will  tend  towards. 
The  area  between  these  curves  represents  the  absolute 
error  calculation  of  R  according  to  eq  2.  The  following 
equation  as  derived  in  Appendix  A  represents  the  per¬ 
centage  of  error: 


lZlrc  (1-e-<y/fc) 

WUZ 


(7) 


where 

Z~q  2-q, 

E%  *  error  percentage  for  a  step  change 
ti  -  time  when  step  input  occurs 
tb  '  time  when  measurement  ends 
q,  *  initial,  steady-state  heat  flux 
q2  s  potential  heat  flux  after  step 
tc  *  characteristic  time  constant  of  building  com¬ 
ponent. 


Figure  9  illustrates  how  the  example  walls  would  in¬ 
fluence  the  measurement  of  R-value  after  a  step  change 
for  a  variety  of  circumstances: 

1.  The  step  in  Afis  up  or  down  10°F. 

2.  The  indoor  temperature  is  always  70°F. 

3.  The  colder  temperature  of  the  step  is  50°  F  or 
-1 0°F. 

4.  There  is  a  variation  in  the  time  after  measure¬ 
ment  began  when  the  step  occurs. 

5.  There  is  a  variation  in  the  time  after  the  step 
when  measurement  stops. 

The  figure  caption  and  key  explicitly  show  all  the  values 
used  in  the  equations  except  that?,  =  (70° -Tw)/R  and 
q2  =  (70 e-Tc)/R  for  a  step  increase  in  AT;  q,  and  q2 
swap  these  values  for  a  step  decrease. 

A  sample  calculation  demonstrates  the  use  of  eq  7 
for  a  2x4  frame  wall.  Assume  f,  =  15  hours,  fb  =  18hr, 
time  constant  =  1 .25  hr,  AT,  =  (70°F-50°F)  and  AT} 

=  (70°F-60°F).  Then  q  =  A T/R  and  R  =  14  imply  that 
q,  =  20/14  and  q2  =  10/14.  Therefore,  substituting 
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Figure  7.  Schematic  rendering  of  corresponding  temper¬ 
ature  and  heat  flow  changes  resulting  from  a  step  change 
in  outdoor  temperature  from  T0  tp  T0 

U  »  1 8-1 5  and  Z  =  0.71  -1 .43  into  eq  7  renders  E% 

=  3.4%  in  agreement  with  the  circled  plot  point  in  the 
top  right  graph  of  Figure  9. 

The  graphs  show  that  the  longer  a  step  input  occurs 
after  accumulating  measured  steady-state  values,  the 
smaller  the  effect  on  the  error.  For  each  wall  it  is 
better  to  stop  measuring  than  to  measure  a  relatively 
short  time  longer.  However,  substantially  prolonging 
the  measurement  ultimately  improves  accuracy. 

Note  that  the  error  increases  with  larger  time  con¬ 
stants.  The  same  measurement  period  results  in  in¬ 
creasing  error  when  the  initial  A T  is  smaller. 

Ramp  change 

A  slow  cooling-off  period  from  dusk  to  dawn  or  the 
passage  of  a  warm  front  may  precipitate  a  slow  change 
in  temperature,  which  has  the  approximate  effect  of  a 
ramp  change  in  AT.  The  qt  will  always  lag  behind  a 
ramp  function  of  qp  as  long  as  it  persists.  Forrester 
(1968)  demonstrates  that  the  lag  shown  in  Figure  10 


will  be  asymptotic  to  a  line  parallel  to  the  ramp  dis¬ 
placed  one  time  constant  after  the  heat  flux  potential 
line. 

The  following  equation,  derived  in  Appendix  C, 
describes  the  approximate  percentage  of  error  for  a 
ramp  change  in  AT: 


F  (1/_2)  (U\Z\-SY) 
*  qytbWZ!  2 

where 


(8) 


U> 
U  = 

z= 
s  = 

Y  = 
£%- 

'a  = 
'b  = 

9 1  = 

'c  = 


h 

V'a 

^2*^1 

u-tc 

s\z\/u 

approximate  error  percentage  for  a  ramp 
change 

time  when  ramp  input  begins 
time  when  measurement  ends 
initial,  steady-state  heat  flux 
potential  heat  flux  when  measurement  ends 
characteristic  time  constant  of  building  com¬ 
ponent. 


Figure  1 1  demonstrates  the  response  to  a  ramp  input 
of  the  same  example  walls  as  before.  It  shows  the  per¬ 
centage  error  in  measurement  according  to  several  as¬ 
sumptions: 

1 .  The  rate  of  change  is  5°  F/hr  up  or  down. 

2.  The  indoor  temperature  is  always  70°F. 

3.  The  colder  temperature  of  the  ramp  input  is  -20°, 
0°,  30°  or  50° F,  depending  on  starting  tempera¬ 
ture  and  duration  of  measurement. 

4.  The  time  between  beginning  of  measurement 
and  ramp  commencement  varies. 

5.  The  time  between  ramp  commencement  and  the 
end  of  measurement  varies. 

We  calculate  and  q2  as  before,  using  70°F,  Tc  and  Tw. 

A  sample  calculation  demonstrates  the  use  of  eq  8 
in  Figure  11  for  a  2x4  frame  wall.  Assume  f,  =  15  hr, 
fb  =  1 7  hr,  time  constant  =  1 .25  hr,  AT  1  =  (70°-50°F) 
and  AT2  =  (70°-60°F).  Then  q  =  A T/R  and  R  =  14 
imply  that?,  =  20/14  and q2  =  10/14.  Therefore, sub¬ 
stituting  U  =  1 8-1 7,  Z  =  0.71  -1 .43,  S  =  2-1 .25  and 
Y  =  (0.75)(0.1 1)12  into  eq  8  renders  E%  =  2.6%  in 
agreement  with  the  circled  plot  point  in  the  top  right 
graph  of  Figure  1 1 . 

The  ramp  change  in  AT  in  Figure  1 1  is  a  more  severe 
case  than  a  step  change  of  the  same  magnitude  in  Figure 
9.  The  top  curves  in  the  key  of  both  figures  represent 
a  jump  in  AT  with  a  change  of  outdoor  temperatures 
between  50°  and  60°F  yet  the  error  magnitude  from 
the  ramp  change  is  consistently  greater. 
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Figure  8.  Step  change  In  AT  and  corresponding  qp.  The  measured  q,.  re¬ 
sponds  In  a  delayed  manner,  creating  the  error  between  the  curves  for  qp 
and  qr . 
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Figure  9.  Percentage  error  from  a  step  change  In  AT  as  a  function  of  time  the  step  occurs  after  measurement 
began.  The  key  shows  the  cold  and  warm  outdoor  temperatures,  Tc  and  Tw,  on  either  side  of  the  step  as  well  as 
the  duration  of  measurement  for  each  wall  after  the  step.  The  wall  properties  are:  2x  4  frame  (R-  1 4),  t cm  t  .25; 
2x6  frame  (R  =  20),  tc  =  1. 78;  masonry  (R  =  7 1),  tc  =  2.82. 
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Figure  10.  Ramp  change  in  AT  and  corresponding  qp.  The  measured  qr  lags  behind  the 
qp  by  one  time  constant,  creating  the  error  between  the  curves  for  qp  and  qr. 


The  graphs  show  that  a  ramp  change  has  a  smaller 
impact  on  measurement  accuracy  die  longer  it  occurs 
after  accumulating  measured  steady-state  values.  For 
each  wall  it  is  better  to  stop  measuring  than  to  measure 
a  relatively  short  time.  However,  substantially  pro¬ 
longing  the  measurement  ultimately  improves  accuracy. 

The  graphs  give  us  a  basic  appreciation  of  the  sig¬ 
nificance  of  changes  in  AT  on  measurement  accuracy. 
Again,  the  larger  tc  Is,  the  greater  the  error. 


APPLICATION  OF  THEORY 

When  we  prepare  to  determine  the  R-value  of  a 
building  component,  we  should  estimate  in  advance 
its  characteristic  thermal  time  constant,  using  eq  3. 
Duration  of  measurement  depends  on  the  temperature 
conditions  on  both  sides  of  the  component  during 
measurement.  Generally,  multiples  of  24  hours  are 
appropriate  durations  for  measurement  of  cyclical 
changes  because  the  accumulated  error  is  likely  to  can¬ 
cel  itself  out.  However,  underlying  warming  or  cooling 
trends  may  impose  themselves  upon  the  diurnal  cycle. 
These  warrant  analysis  as  ramp  changes  according  to 
eq  6  or  Figure  1 1 . 

Overcast  skies  and  precipitation  during  an  approaching 
warm  front  may  dampen  the  effect  of  the  diurnal  cycle 


to  where  the  underlying  warming  trend  is  dominant. 

Such  possibilities  should  prompt  us  to  be  aware  of  the 
weather  conditions  during  measurement  and  monitor 
the  data  to  determine  whether  enough  information  is 
present  for  satisfactory  accuracy. 

In  analyzing  the  data,  we  should  try  to  identify  lags 
between  A T  and  heat  flux  response  qr  to  corroborate 
our  estimate  of  the  building  component’s  time  con¬ 
stant.  Tr,«n  we  should  base  calculations  of  R-value, 
using  eq  2,  on  corresponding  pairs  of  AT  and  qr  where 
the  latter  value  lags  the  first. 

The  application  of  these  theoretical  considerations 
to  actual  field  data  is  imperfect  Figure  1 2  is  a  graph 
of  measured  temperature  changes  and  the  thermal  re¬ 
sponse  of  a  2x4  frame  wall  of  a  building  insulated  with 
urea-formaldehyde  (UF)  foam  recorded  with  the  appa¬ 
ratus  shown  in  Figure  1 .  The  measured  thermal  resistance, 
calculated  according  to  eq  2  using  three  days  of  data, 
is  22  hr  ft2  °F/BTU.  The  theoretical  R-value  from 
published  sources  is  19;  however,  measurement  over 
a  24-hour  cycle  would  have  given  an  R-value  of  25.4, 

19-.9  or  22.5  for  each  respective  day. 

There  is  a  three-hour  lag  between  the  peak  or  valley 
of  A T  over  24  hours  and  dividing  by  the  corresponding 
lagged  heat  flux  does  not  improve  the  results  of  the  in¬ 
dividual  daily  R-values.  The  three-hour  lag  offers  an¬ 
other  difficulty  for  analysis.  According  to  Table  1 ,  a 


11 


Figure  1  1.  Percentage  error  from  a  ramp  change  In  AT  as  a  function  of  time  the  ramp  occurs  after  measurement 
began.  The  key  shows  the  cold  and  warm  outdoor  temperature,  Tc  and  T^,,  at  the  beginning  and  end  of  the  ramp 
as  well  as  the  duration  of  measurement  for  each  wall  after  the  ramp  began.  The  wall  properties  ore  2x4  frame 
(R  *  14),  t,  =  1.25;  2x  6  frame  (R  =  20),  tc  =  1. 78;  masonry  (R  *1 1),  tc  *  282. 
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Figure  12.  Three-day  plot  of  AT  and  heat  flux  on  2x  4  frame  wall. 


Tim*  of  Day 

Figure  13.  Nineteen-hour  plot  of  AT 
and  heat  flux  on  a  2x6  frame  wall. 

fully  insulated  2x4  frame  wall  would  have  a  response 
lag  of  half  that  period.  AT  has  been  converted  to  qp 
using  the  measured  thermal  resistance.  The  qr  shows 
slightly  more  attenuation  than  the  curves  in  Figure  12 
suggest 

Figure  13  depicts  the  thermal  response  of  a  2x6 
frame  wall  containing  UF  foam  and  fiberglass.  The 
expected  R-value  was  24.  However,  the  measured  value, 
according  to  eq  2,  was  only  R  =  16  with  19  hr  of  data. 

If  we  log  the  last  1 0.5  hours  of  data  as  lagging  1 .8  hr 
to  correspond  with  die  estimated  time  constant  of  the 
wall,  the  R-value  is  closer  to  1 8. 


Laboratory  tests  employing  a  Dynatech  Rapid-k 
thermal  testing  machine  also  highlighted  differences 
between  published  values  for  materials  and  experimental 
results.  Published  values  for  all  parameters  except 
specific  heat  agreed  closely  with  those  measured.  Time 
constants  for  1  in.  of  cork  were  expected  to  be  10  min, 
but  1 8  min  was  the  measured  value.  This  contradicts 
the  results  of  Poppendieck  (1976).  For  natural  soft 
rubber  16  min  was  expected  and  40  min  was  measured. 
These  results  are  more  fully  developed  in  Appendix  F. 

The  time  constraints  on  thermal  measurement  of 
building  R  values  are  only  part  of  the  problem.  These 
examples  of  field  measurements  suggest  that  measure¬ 
ment  techniques  need  improvement  and  published 
data  sources  may  not  apply  to  the  material  at  hand.  At 
least  we  can  now  determine  the  effect  of  the  duration 
of  measurement  on  the  accuracy  of  our  results. 
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APPENDIX  A.  PERCENTAGE  ERROR  FROM  A 
STEP  INPUT 


»t»  Dawn  »t»»  Up 


Error  =  shaded  area 

u  u 

f  top  1r)dt=  f  \z\t'tltcdt 


u 

/ 


-t/tc 


=  -tc  \Z\  J  e 


in.-"'  f 


=  fc  \Z\  1-e_t//'c 

Step  down  signal  area 

^  qpdt  =  qy  tb+UZ ; 


Z  =  q2-q i ,  a  negative  value 

Step  up  signal  area 


dt  =  q}  tb+UZ 


Z  =  q2-d\,  a  positive  value 

%  error  »  error/signal  areas 

,  _felZI  1-e_ty/,c 
*  “  d,rb+t/Z 
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APPENDIX  B.  TIME  CONSTANTS  OF  SAMPLE 
WALLS 

Derivation  of  the  technique  for  determining  time 
constants  of  building  components  is  in  Appendix  G. 

Definition  of  variables: 

an  =  rn  dn  c'„  (1  /diffusivity )  where 

rn  =  thermal  resistivity  (hr  ft  °F/BTU) 
c„  =  specific  heat  (BTU/lb) 
dn  =  density  (Ib/fr ) 

gn  =  (an/<ikr  (conversion  factor)  where 
n  =  layer  number 
k  =  layer  chosen  for  normalization 
xn  =  thickness  of  layer  (ft) 

Wood  frame  wall  (2x4  insulated) 

(all  values  from  ASHRAE  1977) 

n  Materia!  rn  dn  en  an  *n 


1  Asbestos  Shingles  3.36  120  0.20  80.6  0.062S 

2  Plywood  sheathing  14.7  27  0.67  266  0.0625 

3  Fiberglass  insulation  38.4  2.0  0.22  16.9  0.2917 

4  Gypsum  wallboard  10.6  50.0  0.26  1374  0.0625 


Time  constant: 
fc  =  1 .78  hr. 

Masonry  wall  (insulated  cavity,  R  =  10.8) 

n  Material _ rn  dn  cn  an  xn 

1  Face  bricks  1.3  130.0  0.190  32.1  0.333 

2  Concrete  Mock  0.259  144.0  0.156  5.82  0.417 

webs 

3  Concrete  block  1.85  144.0  0.156  414  0.167 

faces 

4  Polystyrene  60.0  34  0.29  60.9  0.167 

Insulation  _ _ _ 


Normalize  to  block  web:  dk  =  o2  =  5.82 
0,  =(32.J/5.82)1/2  =  2.35 

g2  =  1 

=  (41 .6/5.82)  ,/2  =  2.67 
gi  =  (60.9/5.82)1/2  =  3.23 
Time  constant: 

7  =  2.82  hr. 


Normalize  to  fiberglass:  ak  =o3  =  16.9 
0,  =(80.6/16.9) 1/2  =2.18 
g2  =  (266/1 6.9) 1/2  =  3.96 
g3  - 1 

04  =  (137.8/16.9)  ^  *  2.85 
Time  constant  of  wall  (from  eq  3): 


=  l  (2.1 8)(0.0625)+(3.96)(0.0625)+(1 )  (0.291 7) 

+(2.85)(0.0625)] 2  (16.9)/tr2 
fc  =  1 .25  hr. 

Wood  frame  wall  (2x6.  insulated) 


n  Material  an  *n 


1  Asbestos  shingles  80.6  0.0625 

2  Plywood  Owe  thing  266.0  0.0625 

3  Fiberglass  Insulation  16.9  0.458 

4  Gypsum  wallboard  137,8  0.0625 


Normalization  to  fiberglass  is  identical  to  before. 
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APPENDIX  C.  PERCENTAGE  ERROR  FROM  A 
RAMP  INPUT 


Kants  Up  Ramp  Doan 


Tim*  Tima 

Figure  Cl.  Error  from  a  ramp  input  (shaded). 


U  =  rb-fa,  *  =  <72-«i  > 5  =  U~h‘ 

Approximate  error  =  shaded  area 

u 

f  (qp-qr)  dt=  M2  U\Z\-U2SY 

Ramp  up  signal  area 


qvdt  =  qx  t^M2UZ-, 


Z  -  q2~q\t  *  positive  value 

Ramp  down  signal  area 


</p<fr  =  9l  tb+M2UZ; 


Z  =  qj-q\ ,  a  negative  value 


%  error  3  error/signal  areas 


APPENDIX  D.  PERCENTAGE  ERROR  FROM  A 
SINUSOIDAL  INPUT 


Figure  D1.  Error  from  a  sinusoidal  input  (area  between  curves). 

All  times  are  from  </p  (0)  =  c.  Values  for  b  and  7",  come 
from  Figures  1 1  and  1 2. 

qp=a  sin  (ir/1 2)  t+c  qr  =  b  sin(rr/1 2)  (r-rL)+c 

Error  =  area  between  curves 


= — \-a  cos  (w/l  2)  t+b  cos(jr/1 2)  (r-f ,  1  2 
jt  nr, 

= —  j  o(cos(*/1 2)  t ,  -cos(*/l  2) 

-b[:os(ir/12)  (f,-rL)' -cos(jr/12)(r2-rL)]J  =  £ 

Signal  area 

B- ^  qpdt  -c  (r2-f,)-l^  x  (cos(*r/12)  f,-cos(*/12)  f2l 

%  error  •  err  or /signal  areas 
£%  * 
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APPENDIX  E.  PERCENTAGE  ERROR  PROGRAMS  FOR  HEWLETT-PACKARD  HP-25  CALCULATOR 


Tit  Ip  Percentage  error:  step,  ramp  A  sinusoid 
Programmer  S.  Flanders _ 


INSTRUCTIONS 


Step  Input _ 

Store  warm  flux 


Store  cold  flux 


Store  time  constant 


Beginning  time 


End  time 


Percent  error 


Repeat  1  -6  or  4-6 


Ramp  Input 


Repeat  1-6  as  above 


Up  &  down  are  reversed 


Sinusoidal  Input 


Convert  to  radians 


8  |  Beginning  time 


INPUT 

DATA/UNITS 


.  Page _ 1 _ of— 1. 


OUTPUT 

OATA/UNITS 


ENTER 


ENTER 


I - If" 


J 

□L  z  j 


1  ^ 


L 


Percentage  error  from  a  step  input 


Page _ 2 — of — i. 


Switch  to  PRGM  mode,  press  0  (  prom  I  ,  then  key  in  the  program. 


|  DISPLAY 

KEY 

ENTRY 

X 

Y 

Z 

T 

COMMENTS 

UNE 

CODE 

00 

•  s,  *•  ' 

x\  \ 

tb 

Step  up  of  A  T 

01 

73  07 

STO  7 

th 

t. 

02 

41 

© 

-u 

03 

23  00 

STOO 

-u 

04 

24  03 

RCL  3 

-U 

■ 

05 

71 

© 

-u/tc 

mu  i 

06 

15  07 

*1 

07 

01 

1 

1 

«i.  ... 

06 

41 

© 

*1_-1 

09 

24  01 

RCL  1 

<7i 

*i-1 

10 

24  02 

RCL  2 

Qi 

<7l 

*,-1 

_ 

11 

41 

© 

-z 

*1-1 

12 

23  06 

STO  6 

-z 

*1-1 

13 

61 

® 

Z(1-*i) 

14 

24  03 

RCL  3 

h 

Z(1-e,) 

15 

61 

® 

E 

E-Z  rr(1-e_y/<c) 

16 

23  05 

STO  5 

E 

17 

24  01 

RCL  1 

E 

18 

24  07 

RCL  7 

th 

<7 1 

E 

19 

61 

® 

<7i  th 

E 

20 

24  00 

RCL  0 

-U 

0,  th 

E 

21 

24  06 

RCL  6 

-z 

-U 

4i  th 

E 

22 

61 

® 

uz 

<7i  tb 

E 

23 

23  04 

STO  4 

uz 

<7i  tb 

E 

24 

51 

© 

<7i  th  &  UZ 

F 

25 

71 

© 

E% 

=^/<7i  th+UZ 

26 

74 

R/S 

Percent  error 

27 

24  05 

RCL  5 

E 

Step  down  of  AT 

28 

24  02 

RCL  2 

q\ 

E 

29 

24  07 

RCL  7 

th 

<7i 

E 

■ 

30 

61 

® 

q\  ^ 

E 

31 

24  04 

RCL  4 

-uz 

tb 

E 

QEEEflHpflNMHB 

32 

41 

© 

q\  tu+UZ' 

E 

33 

71 

© 

tMJJtWVXl  ifiSSSttSf, 

34 

Percent  error 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

■ 

48 

■  B 

49 

_ _ 1 

REGISTERS 

f!TW 

HnHV 

WRISSM 

wttzmm 

*  temp. 

heat  flux 

■mHutUbh 

constant 

R  .UZ 

a.E 

5 

error 

"r“Z 

67  = 

Qn-Qi  _ 

mmm 

I  time 

24 


! 


I  i 

j 

"■  i 

>  .i 


Tit  Ip  Percentage  error  from  a  ramp  input _ 

Switch  to  PRGM  mode,  press  (T)  ,  then  key  in  the  program. 


Page_l _ of _ i. 


25 


Title  Percentage  error  from  a  sinusoid 
Switch  to  PRGM  mode,  press  |J)  (j 


Page _ j—of-  1 


DISPLAY 

KEY 

LINE 

1  CODE  1 

ENTRY 

,  then  key  in  the  prog  im 


COMMENTS 


Input  times 


STO  6 


24  00  I  RCLO 


14  05 


21 


23  07  I  STO  7 


© 


RCLO 


© 


RCL  1 


STO  4 


RCL  3 


RCL  6 


©  © 


RCLO 


cos 


RCL  3 


RCL  7 


I - — 

[RCLO 

© _ 

COS 

f:©_z 
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|©~ 
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© 
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28  _ 4j 
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31  24  02 

32  _ 6j_ 

33  41 

94  24  07 

35  24  06_ 

36  41 

37  24  05 

38  61 

39  24  04 

40  '  41  ' 

41  _ 71 

42  _ 

43 


COS  a* 

«c. 

tr/1 2 

to+t, 

COS  fir 

_“4 _ 

COS  a-> 

<7  C, 

0Sad 

COS  at 

aCi 

C-C 

aCt 

w/12 

C-C 

oC, 

Co 

aC i 

b 

Co 

a  Ci 

b  Co 

aC i 

_ £_ _ 

_ h. _ 

_ h _ 

h-tl _ 

_ c 

jiteishL 


«^M2Hrz+tJL) _ 

C-C  =  COSot,-COS«d 

C,  =  (12/nHC-C) 

E  —  a  C\-b  Co 


jr _ 

Jh. _ 

E _ 

Jz± _ 

E_ _ 

Slh=b.l 

E  1 
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APPENDIX  F.  EXPERIMENTAL  DETERMINATION 
OF  TIME  CONSTANTS 

Laboratory  tests  of  the  thermal  properties  of  cork 
and  natural  soft  rubber  were  to: 

1 .  Corroborate  the  results  of  Poppendieck  et  al. 
(1976). 

2.  Demonstrate  the  time  constant  for  multiple  layers 
of  different  materials. 

3.  Demonstrate  the  thermal  response  to  a  cyclical 
temperature  input. 

In  summary,  the  time  constants  (in  minutes)  were: 
_ Materia / _ Expected  Constant  Observed  Constant 


1  in.  cork 

9.88 

17.8 

1  in.  rubber 

15.5 

40.0 

1  in.  cork  &  1  in.  rubber 

— 

98.0 

These  results  contradict  those  of  Poppendieck  et  al. 
(1 976)  who  expected  an  1 1 .6-hr  time  constant  and 
obtained  a  10-hr  figure  to  "the  ten  percent  value”  with 
a  step  input.  In  this  report  one  time  constant  occurs 
when  a  thermal  response  has  gone  63%  of  the  way  to 
a  step  input  and  2.3  such  constants  bring  the  response 
within  10%.  Poppendieck  et  al.’s  cork  samples  were 
2  in.  thick,  and  therefore  would  respond  one-fourth 
as  rapidly.  Based  on  the  above  results,  the  sample  of 
Poppendieck  et  al.  (1976)  should  have  only  a  2.73-hr 
time  constant,  according  to  their  definition. 

The  observed  results  were  in  reasonable  agreement 
with  the  method  for  calculating  the  time  constant  of 
multiple  layers  of  different  materials  developed  in 
Appendix  G. 

The  cyclical  test  produced  poor  agreement  between 
expected  and  observed  values  of  amplitude  attenuation 
of  the  input  potential  heat  flux  as  an  actual  heat  flux 
response,  as  shown  in  Figure  1 1 .  There  was  excellent 
agreement  between  expected  and  observed  phase  shift, 
using  Figure  1 2.  The  discrepancy  may  be  due  to  the 
fact  that  the  cyclical  input  was  not  a  sine  curve  as 
treated  in  the  figures. 

The  measurement  apparatus  was  a  Dynatech  Corp. 
Rapid  k  thermal  testing  instrument  with  electrically 
heated  and  water  cooled  plates  controlling  the  AT 
across  the  sample.  Temperature  and  heat  flux  sensors 
registered  on  a  strip  chart  recorder. 

The  test  materials  were  cork  and  soft  rubber.  Both 
were  1  ft  square  and  1  in.  thick. 

The  procedure  for  measuring  the  time  constant  of 
a  sample  was: 

1 .  Allow  the  heat  flux  on  the  cool  side  of  die  sample 
to  reach  equilibrium  at  a  low  AT. 

2.  Set  the  warm  side  to  a  substantially  higher 


temperature,  1 09.6°F. 

3.  Continuously  record  the  heat  flux  response  to 
the  new  equilibrium  temperature  value,  AT  = 
37.8®F. 

4.  Determine  time  for  the  response  to  reach  63%  of 
its  final  value. 

Samples  of  cork,  rubber  and  a  cork-rubber  composite 
underwent  this  procedure. 

The  procedure  for  measuring  the  frequency  response 
of  the  cork  and  rubber  sample  was: 

1.  Allow  the  heat  flux  on  the  cool  side  of  the  sample 
to  reach  equilibrium  at  a  high  AT. 

2.  Change  the  setting  hourly  to  the  low  or  high  al¬ 
ternative  temperature  on  the  warm  plate. 

3.  Continuously  record  the  heat  flux  response  to 
temperature. 

4.  Determine  the  attenuation  and  phase  shift  of  the 
heat  flux  response  to  the  heat  flux  potential,  ac¬ 
cording  to  Figures  1 1  and  1 2. 

The  temperature  controllers  on  the  machine  took 
2.2  min  to  achieve  the  criterion  temperature  in  the 
step  temperature  increase.  This  would  cause  a  longer 
response  time  than  a  true  step.  However,  the  AT  over¬ 
shot  the  set  level  by  14%,  which  would  decrease  the 
response  time.  Figure  El  show<  a  plot  from  sample 
data. 

Measurements  of  the  density  and  thermal  resistivity 
of  each  material  agreed  well  with  published  data.  How¬ 
ever,  there  was  no  independent  check  of  the  specific 
heat  of  each  material.  Therefore  the  descrepancy  be¬ 
tween  the  expected  and  the  observed  time  constants 


Figure  FI.  Plot  of  the  heat  flux  response  of 
a  1  -In.  cork  sample  to  a  step  Increase  In  AT. 
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Density  Thickness  Resistivity  Specific  Heat  Constant 
Material  <, Ib/f f)  (ft)  (ft°Fhr/BTU)  (BTU/lb)  {min) 

Cork  14.8*/15+  0.083/NA  32.6/35.7  U/0.48S  17.8/9.9e 

Rubber  62.4/68.6  0.038/NA  12.4/12.5  U/0.48  40/1 5 Je 

e  =  expected,  NA  -  not  applicable,  U  =  no  independent  check,  *  =  Measured 
*  =  Published 


Figure  F2.  Plot  of  the  heat  flux  response  of  combined  1-lh. 
cork  and  rubber  samples  to  a  cyclical  Input  of  tsT  with  a 
2-hr  period. 


remains  unexplained.  The  above  summarizes  the  mea¬ 
sured  and  published  data  for  the  step  inputs. 

The  step  change  determination  of  the  time  constant 
of  rubber  and  cork  together  yielded  tc  =  1 .6  hr. 

The  temperature  controller  took  a  fairly  long  time 
to  cool  down  the  face  plate  with  tapwater  during  the 
cyclical  test.  Consequently,  the  cyclical  inputs  resulted 
in  the  rapid  warming  and  slow  cooling  shown  in  Figure 
F2. 

The  period,  tp  =  2  hr,  of  the  cyclical  input  to  the 
samples  implies  a  time  ratio  of  tjtp  =  0.81 7.  Using 
Figures  1 1  and  1 2,  the  expectations  compare  with  ob¬ 
servations  as  follows: 

Observed 

_ Expected  Value  Value 

Attenuation  [(2irtt/3fp)2+1 |'3/2  *  0.129  1 1/24  =  0.458 

(output/input) 

Phase  lag  (3  arc  un  (2nrcetp)(2  hr/2»)  =  1 .0  1.02  hr 

These  tests  represent  good  indications  of  the  behavior 
sought.  Independent  testing  for  specific  heat  and  fur¬ 
ther  cyclical  testing  with  different  time  ratio  will  help 
explain  unanswered  questions. 


APPENDIX  G.  DERIVATION  OF  TIME  CONSTANT 
FORMULA  FOR  MULTIPLE  LAYERS 

The  exponential  time-dependent  term  in  equation 
4  from  Carslaw  and  jaeger  (1959)  suggest  that  the  time 
constant,  rc  =  ex2/*2,  increases  with  the  square  of  the 
distance  for  a  uniform  slab.  If  we  define  a  -  1  /dif- 
fusivity  and  x  =  thickness  of  slab,  then  a  6-in.  slab  would 
have  a  time  constant  36  times  longer  than  a  1-in.  slab. 

Consider  three  different  slabs  of  the  same  material 
whose  individual  time  constants  are  known  and  com¬ 
pare  their  individual  time  constants  to  die  time  con¬ 
stant  of  the  three  combined  in  one  slab  (Fig.  G1 ). 


o*|  ali^y  ijlV»2 


w2  ir2  w2 

Figure  G 1.  Time  constants  for  Individual  slab 
thicknesses  at  the  same  material  vs  the  time  con¬ 
stant  when  combined  into  one  slab. 

Since  the  chosen  thicknesses  can  be  arbitrary  and  add 
up  to  the  same  total  thickness,  one  constant  D  will 
convert  the  individual  time  constants  to  the  appropriate 
single-slab  figure: 

rc  =  0(o/w2)(x?+x?+x§)  *  (o/«2)(x1+x2+x3)2. 
Therefore, 


*1  *2  *5  *»  *2  ‘J 


Figure  G2.  Time  constants  for  Individual  slab  thick¬ 
ness  of  different  materials  vs  the  time  constant  when 
combined  into  one  slab. 

Choose  a  material,  say  e3,  to  normalize  the  other  pro¬ 
perties  to  according  to  a  constant yn  for  each.  Then, 

a ,  =$i2<»3  fi  *«3  (Pi^i)2/*2 

a2  antl  t2=a3^2x2'>2fn2 

<73=032<73  and  f3  =  o^^x^)2 lit2 ,  where g3  = 

1,  since  this 
material  is 
die  basis  for 
normaliza¬ 
tion. 

Therefore  eq  G1  becomes 

p_fri«rlwy3*e>a  G2 

((Fl*l)2+(02*2)2+fc3Jf3)2l 

In  the  case  of  the  rubber  and  cork  composite  sample, 
xc  =  xr  =  1  /1 2  ft  From  Appendix  F : 


p  _  time  constant  of  single  slab _ 

sum  of  time  constants  of  individual  slabs 

,  (tf/«2)  (x,+x2+x3)2 
W*2)  (x?+x3*x§) 

.  (Xt  »x2+x3)2  (G1 

bi\+x\+x\) 

Suppose  that  each  seynent  is  a  different  material. 
Equation  G1  can  represent  this  if  the  thickness  of  each 
segment  is  adjusted  to  convert  each  separate  material 
into.the  same  equivalent  material  (Fig.  G2). 


a(  *  (0.67  hr )»2/x2  and  ac  =  (0.30  hr)*2/x2 

gr2  =  0.67/0.30  and  gr  =  1 .49.  Therefore  D  for 
the  two  layers  is 

P.[(M9)(1/12H1)C!/12)12 

(1.49/12)2+(1/12)2 

*1.92. 

The  measured  time  constant  of  the  composite  sample 
should  be  1 .92  times  the  sum  of  the  two  constituent 
time  constants.  Therefore,  from  Appendix  F 
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'expected  =  0-92)(0.67+0.30) 

=  1.86  hr. 

This  is  16%  over  the  observed  value  of  1 .6  hours,  in¬ 
dicating  that  the  analysis  is  in  reasonable  agreement 
with  the  observed  value. 

In  conclusion,  the  time  constant  of  a  composite 
material  is: 

'c  =  K/*2)  £  X)  <7n*nj  2  where 
9n  =  («„K)V2  and 

ak  is  the  reciprocal  of  diffusivity  for  a  chosen  material. 
This  is,  of  course,  the  numerator  of  D. 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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